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Abstract. Stump-toed frogs (genus Stumpffia Boettger, 1881 ) are a diverse group of small-bodied frogs endemic to Mada¬ 
gascar. Seven species of this genus occur on Marojejy, a steep massif in northeastern Madagascar. Here we examine the ele¬ 
vational distribution, phylogenetic position, biogeographic origin, and genetic differentiation of this Stumpffia assemblage. 
We show that none of these species is another’s closest relative, but rather they are all independent lineages that probably 
colonised the Marojejy Massif through repeated immigration events. All of the lineages on Marojejy are most closely re¬ 
lated to species south and southwest of the massif, except one lineage, formerly known as Stumpffia sp. Cao 7 , but here as¬ 
signed to S. sorata as a deep conspecific lineage (and referred to as S. cf. sorata), which occurs also in Sorata, 90 km north 
of Marojejy. The species on Marojejy are typically restricted to narrow elevational ranges, but at least two species, S. cf. 
sorata and S. tridactyla , occur over elevations spanning 1000 metres. We assessed the genetic variability of these popula¬ 
tions, and found considerable haplotype separation in fragments of the mitochondrial 16 S rRNA and nuclear Rag-i genes, 
suggesting some disruption of gene flow associated with elevation. We discuss the biogeographic implications of our find¬ 
ings and, based on previously published data, the evolution of non-overlapping bioacoustic parameters among the diverse 
assemblage of Stumpffia species on the Marojejy massif. 

Keywords. Amphibia, Anura, Stumpffia sorata, Stumpffia sp. Cao 7 , Stumpffia tridactyla, candidate species, deep conspe¬ 
cific lineage. 


Introduction 

The Marojejy massif in northern Madagascar is one of 
the hotspots of diversity of stump-toed frogs of the genus 
Stumpffia Boettger, 1881, a Madagascar-endemic group 
of microhylids that contains 40 species after the revi¬ 
sion of Rakotoarison et al. (2017). The largest part of 
the massif is included in Madagascar’s network of pro¬ 
tected areas as Marojejy National Park (Patel 2015). Six 
species of Stumpffia are known to occur in Marojejy Na¬ 
tional Park, all but one apparently being micro-endem¬ 
ic to the massif, and with distinct patterns of elevation¬ 
al distribution: Stumpffia achillei and S. diutissima occur 
only in the lowland forests from 300-750 m a.s.l., while 
S. grandis is apparently only found at around 1330 m a.s.l., 
and S. roseifemoralis also is most common at higher eleva¬ 
tions; on the contrary, S. tridactyla and S. sp. Cao7 [a can¬ 


didate species briefly discussed by Rakotoarison et al. 
(2017) but not taxonomically described], have much larg¬ 
er elevational distributions, with the former occurring 
from 1330-2026 m a.s.l., and the latter from 310-1330 m 
a.s.l. (Rakotoarison et al. 2017). A seventh species of 
Stumpffia occurring at Marojejy, S. sp. Can, is reported 
for the first time herein (see below) and is known from a 
single individual only. 

All of these leaf-litter frogs are typically miniaturised 
and range in adult body size from ca. 9 mm to almost 
30 mm. They usually have small ranges, not only in ele¬ 
vation as summarized above for several species, but also 
in terms of geographic area since most species are known 
from just one or two localities (Rakotoarison et al. 2017) 
- patterns typically found in small-sized amphibians in 
Madagascar (Brown et al. 2016). Their biogeography ap¬ 
parently involved multiple instances of moderately long- 
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distance (>ioo km) dispersal mixed with localised diversi¬ 
fication (Rakotoarison et al. 2017). 

Three main questions arise from the diversity of 
Stumpffia species on this mountain massif: (1) Which 
species are the closest relatives of the massif-endemics, 
(2) does their phylogeny suggest in-situ speciation on the 
mountain itself, and (3) how much gene flow is there in 
species distributed over large elevational ranges on the 
massif? In the present study, we address these questions, 
based on new survey data gathered in 2016 and DNA se¬ 
quences of additional individuals collected during this sur¬ 
vey, to shed further light on the biogeography of frogs on 
the Marojejy massif in particular, and in northern Mada¬ 
gascar in general. 

Materials and methods 

We conducted fieldwork on Marojejy in November-De- 
cember 2016. Field camps and dates were as follows: 
Camp o (14.4463° S, 49.7852° E, 310 m a.s.l.), 14-15 and 
25-27 November; Camp ‘Mantella’ (14.4377° S, 49.7756° 

E, 456 m a.s.l.) 16-17 and 21-24 November; Camp ‘Simpo- 
na (14.4366° S, 49.7434° E, 1325 m a.s.l.) 17-21 November; 
Camp ‘Marojejia’ (14.4350° S, 49.7606° E, 774 m a.s.l.) 30 
November-4 December. Sites above Camp ‘Simpona in¬ 
cluding a site for which we coined the name ‘Buzzard Rock’ 
(14.4406° S, 49.7400° E, ca. 1560 m a.sl.) and the Pandanus 
forest (14.44755° S, 49.7337° E, 2026 m a.s.l.) and areas be¬ 
tween these sites were visited on 19 November 2016. 

Specimens were collected at night or day by searching 
in the leaf litter guided by the calling of males and through 
opportunistic searches. Specimens were anaesthetised 
and subsequently euthanised in MS-222 solution, fixed 
in 90% ethanol and preserved in 70% ethanol. Vouchers 
were deposited in either the Zoologische Staatssammlung 
Miinchen (ZSM) or the amphibian collections of the Men¬ 
tion Zoologie et Biodiversite Animale of the University 
of Antananarivo (UADBA-A). FGZC and ZCMV refer to 

F. Glaw and M. Vences field numbers, respectively. 

Tissue samples were taken by cutting pieces of leg mus¬ 
cle from the euthanised animals and preserved separate¬ 
ly in 99% ethanol. Males and females were distinguished 
based on field observations (calling behaviour) or presence 
of a vocal sac in males, or eggs in females. 

DNA extraction and sequencing 

Genomic DNA was extracted from muscle tissue samples 
preserved in 99% ethanol using a standard salt extraction 
protocol (Bruford et al. 1992). Separate analyses were 
carried out for fragments of two genes, the mitochondri¬ 
al 16S rRNA (16S) gene and the nuclear recombination¬ 
activating gene 1 (Rag-i). Gene fragments were amplified 
via polymerase chain reaction with primers and protocols 
as in Rakotoarison et al. (2015, 2017): 16S was amplified 
with primers 16SL3 (AGCAAAGAHYWWACCTCGTAC- 


CTTTTGCAT) and 16SAH (ATGTTTTTGATAAACAG- 
GCG), with 90 s at 94°C followed by 33 cycles of 45 s at 
94°C, 45 s at 52°C, 90 s at 72°C, and a final extension step 
of 300 s at 72°C. Rag-i was amplified with primers Ragi_ 
Coph_Fi (CGTGATCGGGTAAAAGGTGT) and Ragi_ 
Coph_Ri (TCGATGATCTCTGGAACGTG), with 120 s 
at 94°C followed by 35 cycles of 20 s at 94°C, 50 s at 53°C, 
180 s at 72°C, and a final extension step of 600 s at 72°C. 
PCR products were cleaned with 0.15 units of Shrimp Al¬ 
kaline Phosphatase (SAP) and 1 unit of Exonuclease I (New 
England Biolabs, Frankfurt am Main, Germany) incubat¬ 
ed for 15 min at 37°C followed by 15 min at 8o°C. Purified 
PCR products were sequenced on an automated DNA se¬ 
quencer (Applied Biosystems ABI 3130XL). Sequencing 
reaction (10 pi) contained 0.2 or 0.3 pi of PCR product, 
0.5 pi of BigDye 3.1 (Applied Biosystems, Darmstadt, Ger¬ 
many) and 0.3 pmol of primer. Sequences were checked 
and edited, and heterozygous positions in Rag-i inferred 
in the software CodonCode Aligner 3.7.1 (Codon Code 
Corporation, Dedham, MA, USA). Newly determined se¬ 
quences were submitted to GenBank (accession numbers 
MK452367-MK452379 and MK457424-MK457436). 

Sequence alignment and analysis 

Sequences were aligned with those from previous stud¬ 
ies in MEGA7 (Kumar et al. 2016). For 16S we used the 
dataset from Rakotoarison et al. (2017). We included all 
sequences available for those species occurring in Maro¬ 
jejy, whereas for all other species we included only one se¬ 
quence each. We then calculated a maximum-likelihood 
tree under a GTR+G model as previously determined by 
model testing (Rakotoarison et al. 2017) in MEGA7, with 
2000 full heuristic bootstrap replicates to assess node sup¬ 
port. 

For Rag-i we first trimmed all sequences to equal 
length and removed sequences containing ambiguities 
that could not be interpreted as heterozygotes. We then 
separated sequences into haplotypes using the Phase algo¬ 
rithm (Stephens et al. 2001) as implemented in DNASp 5 
(Librado & Rozas 2009). Phased sequences were subse¬ 
quently used to construct a haplotype network following 
the approach of Salzburger et al. (2011) with the program 
Haplotype Viewer (http://www.cibiv.at/~greg/haploview- 
er) based on a ML tree computed with MEGA7 under the 
Jukes-Cantor model. 

In addition, we also analysed sequences of a different 
stretch of the 16S gene which were obtained by an Illumina 
amplicon approach as described in Vences et al. (2016). 
These short sequences of 201 bp were available for a larger 
number of specimens of S. tridactyla and S. sp. Ca7 from 
different elevations in Marojejy and were thus used to an¬ 
alyse whether populations at different elevations on this 
massif show mitochondrial differentiation despite geo¬ 
graphical proximity. To represent the encountered varia¬ 
tion, consisting of only few mutations, we reconstructed 
haplotype networks as described above for Rag-i. 
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Results 

Diversity and relationships of Marojejy-endemic 
Stumpffia lineages 

A phylogenetic analysis of a DNA sequence alignment of 
598 bp of a fragment of the 16S gene (Fig. 1) confirmed 
the presence of seven species-level lineages of Stumpffia 
on the Marojejy massif (Fig. 2). These include the six lin¬ 
eages already reported by Rakotoarison et al. (2017), 
plus a specimen of Stumpffia sp. Can, newly reported 
here. This specimen was collected by us in 2016 at 1330 m 
a.s.l. on Marojejy. This candidate species was previous¬ 
ly known only from Ambolokopatrika to the southwest 
of Marojejy. While this single and short gene fragment 
is not suitable to reconstruct deep relationships within 
Stumpffia (see the better resolved multi-gene tree in Ra¬ 
kotoarison et al. 2017), it illustrates that the Stumpffia 
from Marojejy are not each other’s closest relatives, as 
each of them (except S. grandis, which has no obvious 
sister lineage) is grouped with high support with other 
species from different sites. 

Three species of Stumpffia on Marojejy are most close¬ 
ly related to species or lineages from Ambodivoangy, a 
low elevation locality ca. 90 km south of the mountain. 
In the following we give their respective genetic diver¬ 
gences as uncorrected pairwise distance in the 16S frag¬ 
ment typically used for DNA barcoding Madagascar’s 
frogs (Vieites et al. 2009) from Rakotoarison et al. 
(2017): S. tridactyla is deeply divergent (10.7% uncorrect¬ 
ed p-distance) from its sister species S. contumelia from 
Ambodivoangy, S. diutissima is fairly strongly divergent 
(8.5%) from its sister, S. pardus; and S. roseifemoralis from 
Marojejy is also strongly divergent (9.6%) from its sis¬ 
ter lineage, an unconfirmed candidate species from Am¬ 
bodivoangy called S. sp. Ca57 by Rakotoarison et al. 
(2017). Stumpffia achillei is sister to S. analanjirofo from 
Nosy Mangabe and the surrounding mainland, ca. 110 km 
south of Marojejy, and these two lineages are separated 
by just 3.0% uncorrected p-distance. The sister lineage of 
Stumpffia sp. Can from Marojejy and Ambolokopatrika 
is S. kibomena from Andasibe, more than 500 km SSW of 
Marojejy, and these lineages are separated by 6.3% uncor¬ 
rected p-distance. This species may therefore warrant de¬ 
scription, but we will deal with its taxonomy and that of 
other red-bellied Stumpffia elsewhere, once that sufficient 
specimens are available for analysis. 

Stumpffia sp. Cao7 is the only lineage whose closest 
relative is found north of the massif, being sister to speci¬ 
mens of S. sorata from the Sorata massif, 90 km NNW 
of Marojejy, separated by 3.3% uncorrected p-distance. 
These lineages belong to a clade (Clade A of Rakotoari¬ 
son et al. 2017) with its centre of endemism in northern 
Madagascar, in contrast to the other species, which all 
belong to clades with their centres of endemism in east¬ 
ern or northeastern Madagascar (Rakotoarison et al. 
2017). 


Taxonomic status of Stumpffia sp. Ca07 

The taxonomic status of Stumpffia sp. Cao7 warrants com¬ 
ment in light of our new sequence data. In our mitochon¬ 
drial tree (Fig. 1), it was recovered as sister to S. sorata 
with high support. The uncorrected p-distance to this lin¬ 
eage of 3.3% in the typically-used barcoding fragment of 
the 16S gene (Vieites et al. 2009) and of 4.2-4.5% in the 
other fragment of the same gene analysed here exceeds 
the threshold of 3% established by Vieites et al. (2009) to 
identify a lineage as a candidate species. A comparison of 
Rag-i sequences (335 bp) revealed a high variation of alle¬ 
les (haplotypes) within S. sp. Cao7 in Marojejy, with up to 
7 mutations among alleles (Fig. 3). Samples from low eleva¬ 
tion on Marojejy had an exclusive allele differing from the 
one found at higher elevation. Haplotypes of S. sorata from 
its type locality (Sorata Massif) differed by only a single 
mutation from the nearest allele of S. sp. Cao7 (from Maro¬ 
jejy Massif, Camp Simpona at 1330 m a.s.l.), but no allele 
sharing was observed. 

Genetic evidence that S. sp. Cao7 is a separate species 
from S. sorata is therefore eqivocal; the two lineages are 
on the one hand clearly distinct in mitochondrial DNA, 
but only at a relatively low level. Given their separation of 
90 km, but connection through more or less continuous 
forest at 1330 m a.s.l., and the absence of surveys from the 
intervening forests, it seems likely that intermediate popu¬ 
lations exist that would connect these two extremes of their 
distribution, and that might break up the 16S distance and 
have allele sharing in nuclear genes. 

We compared the morphology of specimens of S. sp. 
Cao7 to S. sorata, and found them to be extremely sim¬ 
ilar (measurements in Table 1, compared with data from 
Rakotoarison et al. (2017); see also Revised circumscrip¬ 
tion of Stumpffia cf. sorata in Supplementary Appendix). 
The only difference between these lineages that appears to 
be more or less consistent is that specimens from Marojejy 
tend to be slightly smaller than those from Sorata (SVL: 
11.3-15.1 mm vs.15.6-16.0 mm) and have marginally longer 
feet (FOTL/SVL: 0.67-0.76 vs. 0.63-0.65). As bioacoustic 
data are not available from S. sorata, we cannot compare 
the calls of S. sp. Cao7 (described in Rakotoarison et al. 
2017) to those of that species. 

In light of this evidence, to avoid taxonomic inflation, 
we here consider S. sp. Cao7 in a preliminary way as a deep 
conspecific lineage of S. sorata, and refer to it herein as 
S. cf. sorata, until its status can be clarified by further data 
on the genetics, morphology, and bioacoustics of interven¬ 
ing populations between Marojejy and Sorata. 

Elevational distribution and gene flow 
in Stumpffia on Marojejy 

There is considerable variability among species in eleva¬ 
tional range sizes (Fig. 3): S. achillei was found from 480- 
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Stumpffia diutissima ZCMV 15077 Marojejy 
Stumpffia diutissima ZCMV 2082 KC351471 Marojejy 
Stumpffia diutissima ZCMV 2067 KC351470 Marojejy 
Stumpffia diutissima FGZC 2751 EU341075 Marojejy 
Stumpffia diutissima FGZC 2742 KC351469 Marojejy 
Stumpffia diutissima FGZC 2730 KC351468 Marojejy 




■ Stumpffia pardus FGZC 4239 KC351437 Ambodivoangy 

- Stumpffia edmondsi FGZC 2677 EU341076 Vohidrazana 

- Stumpffia tetradactyla ZCMV 3396 KC351405 Nosy Boraha 


■ Stumpffia garraffoi DRV 5621 KC351461 Mahasoa forest 

- Stumpffia spandei ZCMV 8802 KC351472 Mahasoa forest 
— Stumpffia jeannoeli ACP 2251 Betampona 
— Stumpffia nigrorubra FGMV 2002163 KC351411 Ranomafana 

- Stumpffia makira ZCMV 11257 KC351460 Makira Fotsialanana source 
— Stumpffia dolchi ZCMV 885 KC351484 Nosy Mangabe 

Stumpffia sp. Call ZCMV15214 
Stumpffia sp. Call MRSN A2583 AY594082 Ambolokopatrika 
Stumpffia kibomena AND 020 KM881485 Andasibe 

_I Stumpffia grandis FGZC 2820 EU341078 Marojejy 

99 1- Stumpffia grandis FGZC 2828 Marojejy 
— Stumpffia fusca FGZC 4253 KC351369 Ambodivoangy 

Stumpffia analanjorofo ZCMV 2178 KC351365 Nosy Mangabe 
Stumpffia achillei FGZC 2729 EU341079 Marojejy Camp Mantella 
| Stumpffia achillei ZCMV 2095 KC351401 Marojejy below Camp Mantella 
Stumpffia achillei ZCMV15143 Marojejy 
Stumpffia achillei FGZC 2750 Marojejy Camp Mantella 



99 


54 Stumpffia achillei FGZC 2765 KC351400 Marojejy Camp Mantella 
Stumpffia achillei ZCMV 861 KC351402 Marojejy Camp Mantella 
Stumpffia achillei ZCMV 866 KC351403 Marojejy Camp Marojejia 

- Stumpffia betampona ACP 1870 Betampona 

- Stumpffia miery ZCMV 5868 KC351448 Ambohitsara 

■ Stumpffia miovaova FGZC 3675 Sorata 

- Stumpffia obscoena FGZC 2666 EU341073 Andasibe 

— Stumpffia davidattenboroughi ACP 2004 Betampona 

— Stumpffia meikeae DRV 6169 KC351388 Tsaratanana 



99 r 


— Stumpffia hara FGZC 1293 GU989210 Nosy Hara 

- Stumpffia megsoni FGZC 1951 GU989213 Ampombofofo 

■ Stumpffia be FGZC 1643 GU989228 Ankarana 
- Stumpffia Staffordi FGZC 1601 Montagne des Francais 


- Stumpffia analamaina ZCMV 11432 KC351446 27 km NE Antsohihy 

- Stumpffia yanniki DRV 6412 KC351453 southwest of Bealanana 

- Stumpffia madagascariensis FGZC 390 GU989234 Mt Ambre 

— Stumpffia pygmaea MV 2001F21 EU341065 Nosy Be 
■ Stumpffia psologlossa FGMV 2001G42 KC351418 Nosy Be 

- Stumpffia larinki FGZC 1620 GU989206 Ankarana 

971 — Stumpffia angeluci FGZC 436 GU989216 Montagne dAmbre 
56fl—. 

99 


75 r 


i_ 


61 


94 


99 


Stumpffia maledicta ACZC 1644 KC351452 Montagne dAmbre 
Stumpffia huweiZCNN 13620 Montagne dAmbre 
— Stumpffia mamitika FGZC 2725 EU341074 Andrakata-Andapa 

- Stumpffia gimmeli FGMV 2001 .G6 KC351412 Benavony 

— Stumpffia iharana AMNH Al67359 DQ283411 Vohemar - Bezavona 
- Stumpffia sorata FGZC 3617 Sorata 
I Stumpffia cf. sorata (Ca7) ZCMV 15070 Marojejy Camp 0 
1 Stumpffia cf. sorata (Ca7) FGZC 2728 KC351477 Marojejy Camp Mantella 
Stumpffia cf. sorata (Ca7) ZCMV 15050 Marojejy Camp 0 
Stumpffia cf. sorata (Ca7) ZCMV15167 Marojejy Camp Simpona 
Stumpffia cf. sorata (Ca7) FGZC 2826 EU341083 Marojejy Camp Simpona 
9* Stumpffia cf. sorata (Ca7) FGZC 2827 KC351479 Marojejy Camp Simpona 
Stumpffia cf. sorata (Ca7) FGZC 2870 KC351478 Marojejy Camp Simpona 
991 Stumpffia cf. roseifemoralis FGZC 4238 KC351380 Ambodivoangy 
Stumpffia cf. roseifemoralis FGZC 4250 KC351379 Ambodivoangy 
Stumpffia roseifemoralis FGZC 2858 EU341072 Marojejy 
Stumpffia roseifemoralis ZCMV 15172 Marojejy 

_ Stumpffia roseifemoralis FGZC 2047 KC351378 Marojejy 

Stumpffia roseifemoralis FGZC 2808 KC351376 Marojejy 
Stumpffia roseifemoralis FGZC 2883 KC351377 Marojejy 
Stumpffia roseifemoralis ZCMV 2047 Marojejy 

- Stumpffia contumelia FGZC 4248 KC351372 Ambodivoangy 

Stumpffia tridactyla ZCMV 15285 Marojejy Camp Simpona 
Stumpffia tridactyla ZCMV 15290 Marojejy 
Stumpffia tridactyla ZCMV 15284 Marojejy Camp Simpona 
Stumpffia tridactyla FGZC 2844 EU341080 Marojejy 
| Stumpffia tridactyla ZCMV 15227 Marojejy high elevation 
Stumpffia tridactyla ZCMV 15282 Marojejy high elevation 
Stumpffia tridactyla ZCMV 15283 Marojejy high elevation 
- Anilany helenae FGZC 2101 Ambohitantely EU341070 


99 



99 

83 

85 


% 


Figure 1. Maximum Likelihood tree inferred from 598 bp of the mitochondrial 16S gene, with all available samples of Stumpffia occur¬ 
ring in Marojejy (coloured), plus one sequence each of all other nominal species of the genus. Numbers at nodes are support values in 
percent from a bootstrap analysis (2000 replicates; only values >50% are shown). Specimens assigned to S. sorata are shown in a grey box. 
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750 m (restricted to bamboo forests), and S. diutissima 
from 310-750 m. Stumpjfia roseifemoralis was found mostly 
at 1330 m around Camp Simpona, but one record (in need 
of confirmation; see Discussion) exists from at 481 m, and 
is therefore possibly distributed over this wide elevational 
range. Stumpffia grandis and S. sp. Can were found only 
at 1330 m, S. sorata was found from 310-1330 m, and S. tri- 
dactyla from 1330 m to the peak at nearly 2300 m. 

We analysed genetic differentiation related to elevation in 
two of these species, S. tridactyla and S. sorata, the elevation¬ 
al ranges of which both span over 1000 m a.s.l. In S. sorata, 
specimens from low elevations shared a unique haplotype 
that differed by at least two steps from high-elevation speci¬ 
mens in Rag-i sequences (Fig. 4). Additionally, four speci¬ 
mens from 1326 m a.s.l. differed by two mutations from 
three specimens collected between 310 and 481 m a.s.l. in a 
short segment of the mitochondrial 16S rRNA gene (Fig. 5). 
In S. tridactyla, three specimens found at 2026 m a.s.l. dif¬ 
fered by 1-2 mutations from four specimens collected be¬ 
tween 1326-1573 m a.s.l. in the same 16S fragment (Fig. 5). 


Discussion 

Phylogeographic origin of Stumpffia species on Marojejy 

A previous phylogenetic study included various 
Stumpffia species from Marojejy (Wollenberg et al. 
2008) and placed two of these Marojejy lineages into a 
clade (S. grandis and S. sp. 4, the latter corresponding to 
S. achillei), suggesting that these might represent an event 
of in-situ diversification on the massif. The wider analysis 
of Rakotoarison et al. (2017) confirmed by additional 
data herein (Fig. 1) uncovered and described an unprec¬ 
edented number of additional species of this genus, and 
revealed that the sister lineages of the Marojejy Stumpffia 
species all are allopatrically distributed, not occurring on 
this massif. 

Many of these occur to the south of Marojejy, i.e., in the 
Makira-Masoala lowland forest (localities Ambodivoangy 
and Nosy Mangabe) or the adjacent Ambolokopatrika for¬ 
est, while S. sorata is found in Sorata to the north. All of 
these sister lineages occur in the north-eastern or north- 



Figure 2. Stumpffia species occurring in Marojejy National Park: (a) Stumpffia sp. Ca7 (here considered as S. cf. sorata: ZSM 555/2016 
(ZCMV 15182); (b) Stumpffia roseifemoralis, ZSM 529/2016 (ZCMV 15172); (c) Stumpffia grandis (not assigned to voucher speci¬ 
men); (d) Stumpffia diutissima, ZSM 547/2016 (ZCMV 15117); (e) Stutnpffia sp. Call, ZSM 551/2016 (ZCMV 15214); (f) Stumpffia 
cf. sorata, ZSM 544/2016 (ZCMV 15181); (g) Stumpffia tridactyla (specimen collected, but specimen number uncertain); (h) Stumpffia 
achillei, ZSM 536/2016 (ZCMV 15149). 
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Table 1. Morphometric measurements of Stumpffia cf. sorata from Marojejy (all in mm). ND: Not Determined 


Collection number (Field nr) 

Sex 

SVL 

HW 

HL 

TD 

ED 

END 

NSD 

NND FORL HAL 

HIL 

FOTL FOL 

TIBL 

ZSM 379/2005 

(FGZC 2826) 

ND 

14.6 

5.1 

4.9 

1.1 

1.6 

1.3 

0.9 

2.1 

7.9 

2.7 

21.1 

10.0 

6.0 

7.4 

ZSM 537/2016 

(ZCMV 

15047) 

ND 

11.4 

4.2 

4.6 

1.2 

2.0 

1.1 

0.9 

1.4 

7.6 

2.4 

16.8 

8.7 

5.0 

6.4 

ZSM 538/2016 

(ZCMV 

15048) 

ND 

11.3 

4.1 

4.6 

1.1 

1.9 

1.0 

0.9 

1.2 

7.7 

2.3 

16.9 

8.7 

5.2 

6.3 

ZSM 539/2016 

(ZCMV 

15049) 

ND 

14.0 

4.7 

4.8 

1.3 

2.0 

1.1 

1.1 

1.2 

8.5 

2.7 

20.5 

8.9 

6.0 

6.4 

ZSM 540/2016 

(ZCMV 

15050) 

ND 

12.6 

4.8 

4.4 

1.2 

2.0 

1.1 

0.9 

1.6 

6.8 

2.8 

19.3 

8.9 

5.2 

6.3 

ZSM 541/2016 

(ZCMV 

15069) 

ND 

11.6 

4.2 

4.7 

1.2 

2.0 

1.1 

0.8 

1.1 

7.5 

2.9 

16.7 

8.8 

5.2 

6.4 

ZSM 542/2016 

(ZCMV 

15070) 

ND 

12.0 

4.5 

4.7 

1.2 

2.0 

1.1 

1.0 

1.4 

7.7 

2.7 

20.8 

8.7 

5.2 

6.4 

ZSM 543/2016 

(ZCMV 

15168) 

ND 

14.8 

6.0 

4.9 

1.2 

2.2 

1.4 

0.9 

1.9 

6.7 

2.9 

20.7 

9.8 

6.7 

6.6 

ZSM 544/2016 

(ZCMV 

15181) 

Male 

14.6 

5.3 

5.0 

1.2 

2.2 

1.4 

0.9 

1.5 

8.6 

3.2 

22.4 

11.6 

6.8 

7.7 

ZSM 545/2016 

(ZCMV 

15182) 

Male 

15.1 

5.0 

5.2 

1 

2.5 

1.2 

1.1 

1.6 

6.1 

2.9 

22.3 

10.3 

5.8 

6.8 

ZSM 550/2016 

(ZCMV 

15212) 

ND 

15.1 

5.1 

5.2 

1.1 

2.4 

1.2 

1.1 

1.5 

6.0 

2.8 

22.1 

10.1 

5.5 

6.5 

ZSM 554/2016 

(ZCMV 

15281) 

Juvenile 

8.5 

3.2 

3.4 

0.3 

1.5 

0.9 

0.8 

1.1 

4.3 

1.9 

12.7 

5.7 

2.9 

4.0 

ZSM 555/2016 

(ZCMV 

15187) 

ND 

14.4 

5.1 

4.6 

0.9 

2.5 

1.5 

0.9 

1.4 

7.5 

3.0 

ND 

ND 

ND 

ND 

ZSM 556/2016 

(ZCMV 

15167) 

Juvenile 

9.2 

3.4 

3.5 

0.3 

1.7 

0.8 

0.8 

1.1 

4.5 

2.0 

13.9 

6.1 

3.6 

4.4 


ern geographic regions of Madagascar (Brown et al. 2016) 
at distances <200 km, suggesting a mechanism of regional 
diversification (Wollenberg et al. 2008) in allopatry, pos¬ 
sibly in part via divergence in montane refugia (Raxwor- 
thy & Nussbaum 1995). 

Of particular interest is the occurrence of S. sorata both 
in Marojejy (here reported as S. cf. sorata) and the Sorata 
Massif revealed here. We have recently published several 
species descriptions arising from a herpetological survey 


of the Sorata massif, undertaken in 2012 (Scherz et al. 
2015, 2017, 2018a, b, Rakotoarison et al. 2017, Protzel et 
al. 2018). In the majority of these cases, a close affinity has 
been suggested between the taxa of Sorata and Marojejy; 
several species have their sisters lineages on Marojejy (e.g. 
Gephyromantis (Duboimantis) grosjeani sister to G. (D.) 
tandroka; Rhombophryne longicrus sister to R. minuta), 
while others are apparently conspecffic across the massifs 
(e.g. Calumma uetzi, Gephyromantis (Vatomantis) lomo- 



Figure 3. Schematic graph showing known elevational distribution of species of Stumpffia on the Marojejy Massif, along with photos 
of the habitat at various elevations and the elevational profile of the track up the Marojejy summit, with the collection sites mentioned 
in the text and Supplementary materials. 
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rina, G. (Duboimantis) schilfi, and possibly also Rhombo- 
phryne vaventy). The variability in genetic differentiation 
between lineages found on Marojejy and Sorata, from con- 
specific populations with low degrees of genetic differ¬ 
entiation, to distinct species with high degrees of genetic 
and phenotypic differentiation, suggests a long and possi¬ 
bly on-going history of recurrent vicariance and biotic re¬ 
connection of these rainforests. A detailed dating of these 
events is not possible due to the lack of thorough time-cal¬ 
ibrated phylogenies for the majority of the taxa concerned, 
but will be a promising endeavour for future studies. 



Stumpffia sorata 


Figure 4. Haplotype network based on Rag-1 sequences (335 bp) 
of 12 specimens of S. cf. sorata from Marojejy, and 4 specimens 
of S. sorata from the Sorata Massif, showing separation of their 
haplotypes by a single step, and high variation within S. cf. sorata. 
Note that low-elevation specimens of S. cf. sorata (Camp O/Man- 
tella) share an exclusive haplotype. 

S. tridactyla (16S rRNA) 

Marojejy Marojejy 



N =3 N =3 


S. cf. sorata (16S rRNA) 

Marojejy 
Camp Simpona 



Figure 5. Flaplotype networks based on 221 bp of the mitochon¬ 
drial 16S gene. The networks show a low but consistent differ¬ 
entiation between samples of S. tridactyla and S. cf. sorata from 
different elevations on the Marojejy massif, respectively. Note that 
this is a different 16S fragment than used in the phylogenetic tree 
in Figure 1. 


Stumpffia sorata is among the species apparently shared 
between the Sorata massif and the Marojejy massif. Par¬ 
ticularly noteworthy in this case is that S. sorata is not re¬ 
stricted just to high elevations like other species shared be¬ 
tween these two sites, but instead occurs across almost l 
km elevational range in Marojejy. Whether or not the spe¬ 
cies formerly occurred over such a broad range in Sorata is 
unknown, but forest below 1300 m is practically eradicated 
on that massif, so low-elevation populations are probably 
extinct there. They may, however, persist on the neigh¬ 
bouring Andravory massif. Although S. sorata occurs from 
300-1330 m a.s.l. on Marojejy, the lineage genetically clos¬ 
est to that in Sorata is found at 1330 m a.s.l. At elevations 
around 1120 m a.s.l. there is continuous (and forested) con¬ 
nection between these massifs, which constitutes roughly 
the minimum possible single-elevation distance between 
the localities. At higher elevations there is no continuity, 
while at lower elevations, there is continuity, but the dis¬ 
tance travelled is greater, encompassing the perimeter of 
additional, lower massifs. Anthropogenic deforestation has 
also recently broken continuity in some areas, especially 
at lower elevations. Thus, gene flow between populations 
on these two massifs, as in S. sorata, is probably currently 
greatest or at least easiest at 1100-1330 m a.s.l. This pattern 
could be explicitly tested in the future based on population 
genomic data sets and resistance mapping and least-cost 
path analysis (Chan et al. 2011). Major climate transitions 
may have pushed species toward and away from the ‘opti¬ 
mal’ elevations for gene flow, and the timing of divergence 
between sister-lineages from Sorata and Marojejy should 
be compared to identify whether the degree and timing of 
differentiation are correlated. 

Stumpffia sorata is also the only Stumpffia species 
present on the Marojejy massif that belongs to a clade that 
has its biogeographic centre of diversity and probable ori¬ 
gin in northern Madagascar (Rakotoarison et al. 2017). 
Marojejy represents a southeastern expansion of this clade, 
and we consider it likely that S. sorata originated in Sorata 
and arrived in Marojejy via dispersal along the connecting 
mountain chain. 


Bioacoustic and morphological divergence 
in syntopic Stumpffia on Marojejy 

Biogeographically, the available data suggest that the com¬ 
munity of Stumpffia species on Marojejy is due to commu¬ 
nity assembly through independent colonisations rather 
than in-situ diversification. However, the species making 
up this community show a clear differentiation in mor¬ 
phology and bioacoustics, suggesting community assem¬ 
bly may not have been random. Marojejy National Park 
is home of six nominal species (S. tridactyla, S. grandis, 
S. roseifemoralis, S. diutissima, S. achillei and S. sorata ) and 
one candidate species (S. sp. Can). Morphologically, the 
Marojejy Stumpffia species are remarkable by the presence 
of the miniaturised species S. tridactyla (8.6-10.6 mm) and 
the large sized species S. grandis (19.3-23.7 mm). The four 
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other nominal species S. roseifemoralis (16.2-18.4 mm), 
S. diutissima (13.4-20.0 mm), S. achillei (14.6-19.7 mm) 
and S. sorata (11.3-15.1 mm) have more or less overlapping 
body sizes. A similar pattern is seen in Montagne dAmbre 
National Park in north Madagascar, where there is a single 
highly miniaturised species (S. madagascariensis), several 
species of intermediate size that are elevationally segregat¬ 
ed (S. huwei, S. maledicta, and S. achillei), and two larger 
species recently discovered (unpubl. data). 

Differentiation of advertisement calls of the Stumpffia 
species from Marojejy is remarkable. The species that oc¬ 
cur in strict syntopy differ strongly in call duration: S. di¬ 
utissima (53-56 ms) and S. achillei (36-52 ms) have short 
calls, and S. tridactyla (101-198 ms) has an intermediate 
call duration. S. sorata (290-299 ms) and S. roseifemoralis 
(276-280 ms) have a relatively long call duration (Rako¬ 
toarison et al. 2017). Inter-call interval, although strong¬ 
ly influenced by temperature and motivation of the calling 
male (Kohler et al. 2017), also shows extreme inter-spe¬ 
cies differences, with S. achillei having the shortest of all 
values in the genus (507-582 ms), S. roseifemoralis (2891- 
3304 ms) and S. sorata (2764-3250 ms) having very long 
intervals, and S. tridactyla (969-1121 ms) and S. diutissima 
(1775-2200 ms) intermediate intervals (Rakotoarison et 
al. 2017). Temporal distribution of calling may also differ 
among species, but at present data are still too inadequate 
to give a clear pattern. 

The distinct differentiation among the calls of all species 
in temporal variables means that the calls of these frogs 
can be relatively easily distinguished by the human ear, 
permitting, with some experience, the assessment of the 
presence of the different species without requiring the of¬ 
ten time-consuming capture of the frogs themselves. Syn- 
topic anuran species very rarely have similar and almost 
never identical advertisement calls (Kohler et al. 2017). 
In Marojejy, this pattern is likely caused by local process¬ 
es of character displacement into non-overlapping rang¬ 
es of temporal and spectral variables. However, a number 
of other hypotheses could be advanced on this topic, e.g. 
that only species with such distinct calls have been able 
to successfully integrate into the local community, or that 
bioacoustic distinction between frogs instead results from 
their phylogenetic disparity. Testing among these and oth¬ 
er alternatives will only be possible once bioacoustic data 
for the sister lineages of the Marojejy species become avail¬ 
able; so far, these are almost completely lacking. 

Endemism and the effect of elevation on gene flow 
in Stumpffia 

The current species of Stumpffia on Marojejy do not ap¬ 
pear to have sympatrically split from each other on the 
massif. However, the two species reliably found occurring 
over wide elevational bands, S. tridactyla and S. sorata, do 
show a genetic divergence related to elevation. Stump-toed 
frogs tend to have small ranges in terms of physical area, 
as well as narrow elevational distributions (Rakotoari¬ 


son et al. 2017). In Marojejy, most Stumpffia species are re¬ 
stricted to the massif itself and most also occur over only a 
narrow elevational range, according to current knowledge. 
Exceptions are S. cf. sorata and S. tridactyla, and possibly 
also S. roseifemoralis, which occur over ca. 1000 m eleva¬ 
tion difference. However, for S. roseifemoralis, the occur¬ 
rence at low elevations is currently only based on a single 
record (ZSM 373/2005 = FGZC 2808; Rakotoarison et al. 
2017) and requires confirmation. The two species with veri¬ 
fied exceptionally large elevational ranges, S. cf. sorata and 
S. tridactyla, possess divergent, exclusive haplotypes in a 
fragment of the mitochondrial 16S rRNA gene (Fig. 5) and 
nuclear Rag-i gene (Fig. 4) at different elevations. This sug¬ 
gests that gene flow up and down the slope is impeded to 
some degree. A more comprehensive sampling along mul¬ 
tiple elevational transects would be necessary to under¬ 
stand if, indeed, elevational (and thus adaptive) divergence 
is at play here. Alternative explanations that require testing 
include non-adaptive mechanisms such as isolation-by¬ 
distance, or different micro-refuges in which refugial pop¬ 
ulations became genetically divergent in episodes of unfa¬ 
vourable climate, and subsequently expanded their ranges 
across the massif. 
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Supplementary data 

Revised circumscription of Stumpffia cf. sorata , with a list of 
newly assigned material and an expanded diagnosis. 
Supplementary Figure SI. Stumpffia cf. sorata (ZSM 544/2016) 
adult male specimen from Marojejy. 

Supplementary Figure S2. Stumpffia cf. sorata (ZSM 544/2016) 
adult male specimen from Marojejy, in life. 

Supplementary Figure S3. Colour variation of Stumpffia cf. sorata 
on Marojejy. 
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